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ABSTRACT 

We present near-infrared and optical observations of the field of the Anoma- 
lous X-ray Pulsar IE 2259+58.6 taken with the Keck telescope. We derive a 
subarcsecond Chandra position and tie it to our optical reference frame using 
other stars in the field. We find a very faint source, = 21.7 ± 0.2 mag, with a 
position coincident with the Chandra position. We argue that this is the counter- 
part. In the J, I, and R bands, we derive {2a) limits of 23.8, 25.6 and 26.4 mag, 
respectively. As with 4U 0142+61, for which a counterpart has previously been 
found, our results are inconsistent with models in which the source is powered 
by accretion from a disk, but may be consistent with the magnetar model. 

Subject headings: stars: neutron — X-rays: individual (IE 2259+58.6) — X-rays: 
stars 



1. Introduction 

IE 2259+58.6 is a member of a class of X-ray pulsars whose energy source is not un- 
derstood, called Anomalous X-ray Pulsars (AXPs, see Israel et al. 2002, for a review). 
Their rotational energy loss is insufficient to explain their X-ray luminosities and there is no 
evidence for a companion star from which matter could be accreted. 

In part, the uncertainty about AXPs arises because of a lack of counterparts at other 
wavelengths. Therefore, we have initiated a program of deep optical and near-infrared ob- 
servations of those AXPs which are not too highly absorbed. This program has led to the 
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detection of the first optical counterpart of an AXP, namely the counterpart to 4U 0142+61 
(HuUcman ct al. 2000a). Our earlier observations of IE 2259+58.6 set stringent limits 
{R > 25.7, I > 24.3 mag, HuUeman et al. 2000b). 

In this Letter, we present in §2 deeper optical observations and new near-infrared obser- 
vations of IE 2259+58.6. In §3, we derive an accurate position from Chandra observations 
(see the companion paper by Patel et al. 2001), which we use in §4 to identify a possible 
near-infrared counterpart to IE 2259+58.6. 

2. Observations 

2.1. Optical observations 

Deep R and I band images of the field around IE 2259+58.6 were taken on the night 
of 2000 September 3, with the Low Resolution Imager Spectrograph (LRIS, Oke et al. 1995) 
mounted on Keck I. The night was photometric and the seeing was 0'.'6 (0'.'9) in R (I). Six 
(twelve) images were averaged for a total integration time of 1 hour in R (I). Less deep, 
wider field of view, B, V, R, and I images were taken at the 60- inch telescope on Palomar 
mountain on the night of 2000 November 18 and were used for the astrometry (§3). 

The reduction was done using the ESO-MIDAS software package, in the same way as 
described by HuUeman et al. (2000b, a). The stacked Keck I band image is shown in Fig. 1, 
the Keck R band image is not shown. It is very similar to the image shown in HuUeman 
et al. (2000b). 

We measured instrumental magnitudes on the Keck images using point-spread-function 
(PSF) fitting as implemented in the DAOPHOT package, and cahbrated these using the 
cahbration found earher (HuUeman et al. 2000b). We estimate single-trial 2a detection limits 
of 26.4 (25.6) mag in R (I) by adding 8.1 (7.6) mag to the mean instrumental magnitude of 
objects that have a PSF fitting error of 0.3 mag^. 

2.2. Near infrared photometry 

Near-infrared images in the J and Kg filters were obtained on the nights of 1999 June 
23 and 24, with the Near Infrared Camera (NIRC, Matthews & Soifer 1994), mounted on 



^The detailed simulations described in HuUeman et al. (2000b) convinced us that this procedure gives 
correct estimates. 
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the Keck I telescope. During both nights observing conditions were photometric. During 
the first night, the seeing was mediocre, varying between (/.'G and 1'.'5, but during the second 
night it was excellent, at &!5 in both Kg and J. Therefore, we only use the data taken on the 
second night. These consist of a total of 36 frames in each of J and Kg, with each (J) 
frame composed of 5 (3) coadded exposures of 12 (20) seconds. 

The data were reduced using MIDAS. After subtracting dark frames, a fiatfield and a 
mosaic were created for each filter in the following way. First, each individual frame was 
normalized to unity. Next, at each pixel the variance over all the (normalized) frames was 
computed. For a field which is not too crowded, the median of these variances should be a 
good approximation of the expected variance of a pixel which in each frame only contains sky. 
Then, for each pixel we used an F-test to determine if its variance was significantly higher 
than the median variance, where we set the level of significance at 0.001. If the variance was 
not significantly higher, we simply took the mean of all frames at that pixel. Otherwise, we 
discarded the highest values at each pixel until the variance was consistent with the median 
variance. From the values left, a new mean was calculated. 

After division by this fiatfield, we used the prescription of Liu & Graham (2001) to 
correct each frame for a persistence effect in the read-out amplifiers and subtracted the 
background, which was determined from the frame itself. 

A few well detected stars were used to determine the spatial offsets between individual 
frames. Next, at each sky position the median and standard deviation were computed 
(excluding physical pixels which were obviously bad). If at any sky position an individual 

frame differed from the median by more than three times the standard deviation, it was 
discarded and a new standard deviation was computed. This was repeated, until the values 
of all remaining frames were within three standard deviations of the median. The resulting 
mosaics, shown in Fig. 1, have the corresponding mean value at each sky position. 

We used DAOPHOT to measure instrumental magnitudes on the mosaics, and calibrated 
these using the infrared standard star SJ 9182 (Persson et al. 1998). In Table 1, we list 
positions and magnitudes of selected objects (see Fig. 1). We also estimated single-trial 2a 
detection limits of 23.8 mag in J and 22.4 mag in Kg in the way described above. 

2.3. Optical Astrometry 

We determined coordinates relative to the USNO-A2.0 catalog (Monet 2000). For the 
60-inch images, 1107 USNO-A2.0 stars overlap. We determine centroids for 942 objects 
that were not saturated and appeared stellar, corrected these for instrumental distortion 
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(J. Cohen, 1997, private comm.), and solved for zero-point position, scale, and position 
angle. After rejecting 193 outliers, the rms residuals are 0'.'30 in both right ascension and 
declination. Therefore, the accuracy of the tie to the USNO-A2.0 system is very good, with 
a formal uncertainty of O'.'Ol (standard deviation of the mean). The real uncertainty may be 
limited by uncertainties in the distortion correction; but we believe it should still be better 
than (/.'OS. 

The optical Keck images arc too deep and the infrared images cover too small a field for 
a direct tic to the USNO-A2.0 catalog. Instead, for those, we used stars in common with a 
short R-band exposure, for which Hullcman et al. (2000b) had determined coordinates in the 
USNO-A2.0 system. For all bands, the tie to the short R-band image is accurate to better 
than 0'/02, and the final tie to the USNO-A2.0 system should again be better than O'.'OS. 

3. The X-ray position 

IE 2259+58.6 was observed with Chandra on 2000 January 11, using the Advanced CCD 
Imaging Spectrometer (ACIS). The analysis of this observation is described in the companion 
paper by Patcl ct al. (2001). Here, we use only the data collected in Timed Exposure mode 
to derive a position for IE 2259-1-59.6. We used the data that were "reprocessed" on 2000 
October 31 through the standard pipeline. This reprocessing corrected minor errors in the 
aspect solution compared to the original processing. Further, in a calibration update made 
on 2001 July 26 the ACIS pixel size was changed from 0.492 arcsec to 0.49131 arcsec. We 
applied this correction manually. 

For IE 2259-1-58.6, which is rather bright, the data suffer from pile-up of photons in the 
image core, as a result of which the source looks ring-shaped with a hole in the center. We 
modeled the data using a Gaussian multiplied by a hyperbolic tangent in radius, scaled to 
approach zero at r = 0.0. The resulting best-fit centroids correspond, using the Chandra 
aspect solution, to a J2000 position on the sky of a = 23'^01°'08^26, 6 = +58°52'44'.'9. The 
uncertainty is limited by systematic effects, to a circle with ~ 0.7" radius (Aldcroft et al. 
2000). 

To improve this position, we tried to identify optical counterparts to other sources 
in the field and use these to obtain a boresight correction. For this purpose, we use all 
seven other sources detected by the source-finding algorithm in the central (S3) chip (see 
Table 2). Three of the sources (XI, X5, and X6) have X-ray spectra with a soft component, 
characteristic of stellar sources with little absorption. All of these have counterparts on the 
Digitized Sky Survey (DSS) plates and positions in the USNO-A2.0 catalog. Object XI is 
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RPl of Rho & Petre (1997); van den Berg & Verbunt (2001) found a G7-K1 IV-III star with 
V = 11.8 ± 0.1 coincident with RPl. The expected X-ray luminosity of the star agrees with 
the one measured for RPl. The stars coinciding with object X5 and X6 are relatively bright 
as well (USNOA2.0 blue and red magnitudes of 17.5 and 16.6 (X5) and 17.7 and 15.3 (X6)); 
given the small probability of finding stars this bright at the X-ray position by chance, we 
conclude that these most likely are the optical counterparts. 

Interestingly, no counterparts are seen on the DSS for any of the other sources, which 
have harder X-ray spectra. However, for X2, which has a low absorbing column, we do 
find a faint blue star on our Palomar images. Given the unusual color, we believe this is 
the counterpart. Furthermore, for X7, which appears absorbed, we find a faint extended 
object at the edge of our 1997 R and I-band LRIS images. Given the crowdedness of the 
field, we cannot be sure this is the counterpart; if it is, X7 is most likely a highly reddened 
extragalactic object. Finally, source X4 also overlaps with the 1997 LRIS images, but is in 
a very confused region. 

In Table 2, we fist positions inferred from our Palomar images, from the USNO-A2.0 
catalog, as well as positions from the DSS (both epochs) using astrometric cafibrations 
relative to the USNO-A2.0 catalogue. None of the objects appear to have significant proper 
motion. We used XI and X6 to establish the boresight correction and X2 to verify our 
solution. We did not use X5 because it is faint and relatively far off-axis in the ACIS image, 
and it is not covered by our Palomar data. 

In the Chandra images, XI, X2, and X6 are off-axis and thus appear elongated, which 
may bias the position derived by the source-finding algorithm. For each X-ray source, we 
computed a model PSF that was rotated by the spacecraft roll angle to match the image 
orientation and fit this model to a two dimensional elliptical Gaussian. The offset of the 
center of this Gaussian from the origin tells us the expected size of the image shift introduced 
by the Gaussian fit. We then fit a Gaussian with the same ax, (Jy, and inclination to the actual 
data and subtract the offset determined from the model PSF. We also fist the uncertainty 
in the centroid location estimated from the fit to the actual data. 

Armed with both optical and X-ray coordinates for these objects, we adjusted the 
Chandra coordinate system using the constant offsets, Sx and 5y, that gave the best agree- 
ment with the Palomar positions. The resulting best J2000 position for IE 2259-1-58.6 is 
a = 23*^01^08^295, S = +58°52'44'/45. We stress that this position is tied directly to our 
optical images, so the only source of systematic uncertainty in the position is the uncertainty 
in the offset, of C.'IQ in each coordinate. This dominates over the statistical uncertainty, and 
corresponds to a 99% confidence error radius of 0'.'60. The corrected Chandra position is 
indicated in Fig. 1. 



-6- 



4. Identification of a Counterpart 

We searched for possible counterparts within the Chandra error circle. We did not 
find any counterpart in our R, I and J images and place upper limits of 26.4, 25.6 and 
23.8 mag, respectively. The R and I limits are nearly a magnitude better than our previously 
published values (HuUeman et al. 2000b). However, in the new near-infrared Kg image, we 
find & Kg — 21.7 ± 0.2 mag object (star 1) within the Chandra error circle (Fig. 1). 

Another object, star 2, while within the ROSAT error circle (see HuUeman et al. 2000b) 
is offset from the Chandra position by 0^85 ± 0'.'25, where the error includes the uncertainty 
in the optical position as well as the systematic uncertainty in the X-ray to optical transfor- 
mation. It is inconsistent with being the near-infrared counterpart of IE 2259-1-58.6 at the 
99.9% confidence level. 

It is possible that star 1 is an unrelated background object. The density of such faint 
objects is 0.04arcsec~^. Hence, thanks to the exquisite astrometric precision of Chandra, the 
probability of chance coincidence is small, less than 3%. 

Thus, most likely we have identified the near-infrared counterpart to IE 2259+58.6. 
Patel et al. (2001) report Nh = 9.3±0.3xl02i cm-^ which, using Ay = jj^j^^ mag (Predehl 
& Schmitt 1995), corresponds to the following extinction values: — 4.3, Aj — 3.1, 
Aj — 1.4, Ak — 0.6 mag, which should be accurate to at least a factor of two. Adopting 
these values we show the optical through X-ray photometry of IE 2259-1-58.6 in Fig. 2. Also 
shown are the blackbody and power-law components of the best-fit model to the X-ray data 
(Patel et al. 2001). 

As for 4U 0142-1-61, we find an extremely large X-ray to optical/infrared fiux ratio, 
which is hard to square with most models proposed for AXP. Models in which the source is 
powered by accretion from a fossil disk (e.g., Chatterjee et al. 2000) were already excluded by 
our previous optical limits. However, this was for a standard disk in which optical emission 
is due to reprocessing of the X-ray fiux. 

HuUeman et al. (2000a) showed that for 4U 0142+61 a disk with an inefficiently radiating 
inner part was excluded as well, but that a truncated disk, such as might be present in a very 
tight binary, was still possible. For IE 2259+58.6, we can also exclude this possibility, since 
the K band fiux and optical limits are inconsistent with the predicted roughly Rayleigh- 
Jeans spectral-energy distribution. For the same reason, we can exclude the possibility that 
IE 2259+58.6 is a hot white dwarf. 

The only model that is without problem, although mostly for lack of predictions, is 
that of Thompson & Duncan (1996), in which AXP are magnetars, powered by the decay 
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of a superstrong, ~ 10^^ G, magnetic fields. We hope that our work will stimulate further 
study into these models. In this respect, we are encouraged by the work done on the X-ray 
emission by e.g. Perna et al. (2001) and Ozel et al. (2001) and the qualitative agreement 
between observations and rough estimates of the expected optical flux made by C. Thompson 
(2001, private comm.). 

With two counterparts identified, observational progress should be rapid. It should be 
relatively easy to determine the broad-band spectral energy distributions of both sources 
and identify other AXPs. More difficult, but more rewarding potentially, would be searches 
for optical pulsations and polarisation. 
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Table 1: Positions and magnitudes of selected objects 
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Note. — The errors quoted are the formal fitting errors. The zero point errors are 0.05 mag in both J 
and Kg. 
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Ta.l)le 2: Clutndm X-ray Point Sources on S3 and tlioir (tentative) counterparts 



source^ 


Q!J2000 




r 

<'J2000 






XL (i.y, oiUj 


00 nn 00 
2o UU 00. 


ooy 


1 CO C^O A^ 
+ 00 OZ 4:1. 


.iz 


(J. U9 


1425-i44oDo4o 


33. 


,396 


A ^ 

47, 


.16 


0. Z 


"P\CC 1 d 

USSi 


33. 


,393 


A ^ 

47, 


00 

.23 


0.2 




00. 


,OOU 


47 


r\A 
.u^ 


u. z 


P60 


33. 


,402 


46. 


.96 


0'.'03 


X2 (1.4, 297) 


23 00 43. 


,360 


1 r r r\ on 

-|-58 50 30, 


A A 

.44 


U.06 


roO 


43. 


,373 


OA 

30, 


.11 


0.09 


X3 (I'.'l, 66) 


23 00 46. 


,351 


+58 53 08, 


.76 


O'.'U 


X4 (1'.'2, 32) 


23 00 55. 


,044 


+58 54 59, 


.05 


0'.'25 


X5 (3'.'5, 37) 


23 00 59. 


,943 


+58 59 11, 


.03 


0'.'72 


1425-14449773" 


59. 


,798 


11, 


.68 


0'.'2 


DSSl ^ 


59. 


.816 


11 


.51 


0'.'2 


DSS2 d 


59. 


.815 


11 


.69 


0'.'2 


X6 (2'/3, 380) 


23 01 00. 


.118 


+58 57 31, 


.01 


O'.'U 


1425-14449921" 


00. 


.032 


31 


.12 


0'.'2 


DSSl ^ 


00. 


.038 


31 


.00 


0'.'2 


DSS2 <i 


00. 


,107 


31 


.20 


0'.'2 


P60 


00. 


,105 


31, 


.17 


0'.'03 


X7%a!9, 27) 


23 01 07. 


.175 


+58 49 41, 


.62 


0'.'09 



"For the X-ray sources we give in parentheses we give the approximate image size and number of counts 
detected. 

''XI and X2 arc RPl and RP2 of Rho & Petre (1997) respectively. 
'=USNO-A2.0 identifier. 

''Epochs are 1953.830 for DSSl and USNO-A2.0, 1989.668 for DSS2 and 2000.885 for P60. 
'^There is an extended source at this position on the edge of a plate taken with LRIS in 1997. 

Note. — All sources except X5 were on the Palomar images. Except for XI, X2 and X6, no object is 
found at the position of any of these sources on these images. All X-ray positions were shifted by +0'.'28 in 
RA and — 0'.'45 in DEC to match the optical positions. The optical positions arc on the system defined by 
the USNO-A2.0 catalog. Hence, they should be on the International Celestial Reference System to ^ 0'.'2. 
Both the DSSl and the USNO-A2.0 positions are determined from the same plate. See text for details. 
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Fig. 1. — Cut outs of the I band image and Kg and J band mosaics of the field around 
IE 2259+58.6. A segment of the Chandra 99% confidence error circle is overplotted on the 
Kg band image, on all other images the position is indicated by tickmarks. The best ROSAT 
error circle of Hulleman et al. (2000b) is overdrawn on the 1 band image for reference. We 
note that although star 1 appears to consist of two dots in Kg, it is consistent with being a 
single point source. 



Fig. 2. — Broad band emission spectrum of IE 2259+58.6. Denoted by CXO are the X-ray 
data (Patel et al. 2001), plusses are absorbed and diamonds unabsorbed X-ray fluxes. Shown 
as dashed lines are the blackbody and power-law components of the best fit model. Also 
shown are the optical and infrared limits and the Kg band detection, as observed and after 
correction for interstellar reddening. 
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